Extended interface states enhance valley splitting in Si/SiC>2 
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Interface disorder and its effect on the valley degeneracy of the conduction band edge remains 
among the greatest theoretical challenges for understanding the operation of spin qubits in silicon. 
Here, we investigate a counterintuitive effect occurring at Si/SiC>2 interfaces. By applying tight 
binding methods, we show that intrinsic interface states can hybridize with conventional valley 
states, leading to a large ground state energy gap. The effects of hybridization have not previously 
been explored in details for valley splitting. We find that valley splitting is enhanced in the presence 
of disordered chemical bonds, in agreement with recent experiments. 
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Introduction. — The transistor revolution has granted 
silicon heterostructures a special status amongst materi- 
als platforms. Yet after many years of intense study, this 
system still reveals new and intriguing features. This 
is due in part to technological advances that open the 
door to new physical regimes. However, the interest in 
Si has also been stirred by its unusual materials proper- 
ties. The Si conduction band (CB) possesses six degener- 
ate minima in the first Brillouin zone, known as valleys. 
Quantum well confinement and/or the application of uni- 
axial strain (e.g., in the case of Si/SiGe heterostructures) 
in the [001] direction reduces the bulk, cubic symmetry, 
and raises the energy levels associated with the transverse 
x and y valleys flf. At very low temperatures, the CB 
physics is therefore governed by the spin and the z valley 
degrees of freedom. Control over valley degeneracy is a 
key concern for Si spin qubits 

Experimental 0, @] and theoretical 0, 0, H| investi- 
gations of the physical mechanisms of valley coupling re- 
veal that a sharp interface between a quantum well and a 
quantum barrier (most commonly SiGe or Si02) can pro- 
duce a sizable energy splitting between the valley states, 
and that roughness can suppress this effect 0, [l(| ■ Real- 
istic theoretical estimates for the interface-induced valley 
splitting are on the order of 0.1-1 meV [ll|, in agreement 
with many experiments. However, they cannot explain 
the recent puzzling results of Takashina et al. 5] . In an 
asymmetrically grown Si/Si02 quantum well, they ob- 
serve a large ground state gap of 23 meV at the buried 
oxide (BOX) barrier, but a more typical valley splitting 
at the second, thermally grown oxide barrier [121 ] . 

In this work, we demonstrate that conventional CB 
electron states tend to hybridize with intrinsic Si/Si02 
interface states (IS), which form in the gap. Hybridiza- 
tion can produce a conducting ground state that is non- 
degenerate, due to strong valley orbit coupling. The re- 
sulting ground state gap can be tens of meV larger than 
the valley splitting between pure CB states, which could 
explain the large values measured in Ref. 5. Such a gap 
would be ideal for controlling spin qubits in Si. 



Tight Binding Model. — Since IS are linked to the 
atomic details of the interface, their theoretical descrip- 
tion should likewise be atomistic. Here we study the 
formation of IS within a minimal, single-electron, two- 
band tight binding (TB) model [llj designed to capture 
the low-energy physics of the CB. We consider both one- 
dimensional (ID) and two-dimensional (2D) versions of 
the model. In the ID version, we consider only hopping 
along the z direction. We extend the ID model of Ref. [ill 
to describe the Si02 barrier as an effective linear chain, 
as sketched in Fig. 0Ja). The full chain is composed of 
three regions: region I, with parameters corresponding 
to bulk Si; region II, involving two sites on either side 
of the interface; and region III, with parameters mim- 
icking the Si02 barrier. For the 2D model, we include 
additional, lateral hopping terms, arranged on a square 
lattice [l9j |. To probe the interface physics, we solve the 
full TB hamiltonian and we analyze the resulting low en- 
ergy spectrum and eigenstates. Details about the models 
and solutions are presented below. 

In region I, the nearest and next-nearest neighbor 
vertical (z direction) hopping parameters (ui and v\) 
are chosen to reproduce the essential features of the 
bottom of the CB: the Si longitudinal effective mass, 
mi = 0.916too, and the CB minima, fco = ±0.82(27r/asi), 
where ast is the length of the Si cubic unit cell. This gives 
u\ = usi = 0.68 eV and vi = vsi = 0.61 eV. For the ID 
model, we adopt the onsite energy t\ = esi = 1-41 eV, 
so that the valley minima occur at zero energy. We also 
add an electrostatic potential, —eFz, to the onsite pa- 
rameter, where typical, experimental electric fields fall in 
the range 0.01 < F < 0.1 V/nm. For the 2D model, we 
introduce an additional nearest neighbor lateral hopping 
parameter, u V j = = —10.91 eV, which gives the cor- 
rect transverse effective mass rrit — 0.191 too- We also 
modify the onsite parameter e\ = e^f — 23.23 eV to cor- 
rectly set the valley minima to zero. 

In an empirical TB model, we determine the hopping 
parameters by matching bulk materials properties. The 
materials properties of the interface are not well known, 



2 



Barrier 





9.0 




6.0 




3.0 


> 


0.0 








>. 




o> 




ci3 


9.0 


c 




LU 


6.0 




3.0 




0.0 




(c) 










L Interface 





Proj. LDOS (eV 



0.05 
0.04 
0.03 
0.02 
0.01 




-2-10 1 
z(nm) 

FIG. 1: (Color online) (a) TB model parameters along z for 
a Si/SiC>2 interface, (b) The projected local density of states 
(pLDOS), nP ro \ as a function of vertical position and energy. 
The interface hopping parameters are taken to be u\\ = usi 
and v\\ = t>si- (c) The pLDOS obtained for hopping parame- 
ters mii = usio 2 an d fn = «Si0 2 • An IS appears in the gap, as 
indicated by a dip in the CB. For panels (b) and (c), the re- 
sults are obtained for a 2D system where model (ii) is adopted 
in region III. 



however. At a fundamental level, the atomic orbitals are 
locally modified by the inhomogeneous chemical environ- 
ment [jjjj . Such variability is amplified at a Si/SiC>2 in- 
terface, where important charge transfer occurs alon g th e 
Si-0 bonds [3] and the interface is highly disordered [3] . 
Thus, while region II is key to understanding IS, the TB 
hopping parameters are not known with certainty. To ad- 
dress this problem, we consider a thorough, continuous 
range of parameters u\\ and vu, as described below. 

Although region III is better understood than region II, 
TB models of SiC>2 are still not well established. We 
therefore explore three minimal models for region III, 
to observe how new behaviors emerge, not aiming for 
a quantitative description. Our results do not depend 
qualitatively on the particular oxide model, as will be 
shown. The three models we consider are given as follows. 

(i) An indirect gap crystalline material. In this case, we 
use the same hopping parameters as Si, so that um = ui, 
I'm = «i, and u\ n = u\. Only the band edge is shifted 
upward, with the experimental offset of W rj 3.0 eV. The 
onsite parameter is then given by em = cj + W. 

(ii) A direct gap crystalline material. In this case, 
there is a single valley at k = 0, and we adopt the 
effective, isotropic mass of 0.34 mo, corresponding to 
FD 3 M /3-cristobalite This is achieved by setting 
Win — «sio 2 — 3.28 eV and vm — 0. The shifted band 
edge is set by em = esio 2 +W — 9.56 eV. 

(iii) An amorphous material. There is strong evidence 
that this is the most accurate model scenario 1141. al- 



though crystalline phases of Si02 have been reported 
within a few layers of the interface with direct 

gap /3-cristobalite being the most abundant phase. Our 
amorphous model consists of TB parameters chosen ran- 
domly in the range uni G [i«si0 2 ~ 0.5, usio 2 + 0.5] and 
win G [—0.5,0.5], in units of eV. For simplicity we take 
the onsite energies and transverse hoppings to be the 
same as model (ii). Thus, the lateral disorder is medi- 
ated entirely by the vertical hopping parameters. 

Interface States. — Intrinsic IS may occur even at 
smooth, ordered interfaces. This is in contrast with ex- 
trinsic IS, which require broken bonds, impurities, or 
some type of disorder potential. The intrinsic IS are of- 
ten referred to as Tamm/Shockley states 
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now demonstrate how they can emerge naturally from 
our Si/Si02 TB models. Later, we will study the depen- 
dence of IS on specific interface models. 

The electron affinity varies abruptly at a sharp in- 
terface, corresponding to a sudden change of the onsite 
terms in our TB model. On the other hand, chemical 
bonds are modified by their local environment, leading 
to a more gradual variation of the hopping terms over 
several monolayers. The energy of the lowest orbital level 
therefore varies locally. This is the origin of the IS [13] • 

To demonstrate the emergence of IS in our TB model, 
we first compute the local density of states (LDOS), de- 
fined as n(i,j,E) = — —Im[G(i, j; i, j; E)], for the co- 
ordinate indices The Green's function matrix 
G(i,j;i',j';E) is obtained from G{E) = l/(E-H + i0). 
In the 2D case, we then compute the projected LDOS 
(pLDOS) along z, defined as n proj (j,£) = T,Mi,j,E). 

The pLDOS provides a means to visualize the local 
CB, as shown in Figs. QJb) and (c). For demonstration 
purposes, we consider two cases. In Fig. [T^b), we take 
Uu = usi and vu = vsi- The CB appears as a dark 
shaded region, with a simple step in the band edge at 
the interface, z = 0. In Fig. [He), we take wtt = wsi0 2 
and vu = Wsi0 2 - I n this case, the lower edge of the CB 
dips by over 200 meV in region II. The dip signifies an 
IS, which is strongly peaked at the interface, and decays 
over a few A on either side. The localized nature of the 
wavefunction suggests that the computed IS is sensitive 
to the values used for uu and vu- 

Ground state gap. — Fig.[2]shows the calculated ground 
state energy gap A, as a function of the interface hopping 
parameters uu and vu- Results are shown for each barrier 
model, (i)-(iii), in region III. In each case, we note that A 
is determined principally by the hopping parameter uu, 
rather than vu- For small un, A ~ 7 meV corresponds 
to the conventional valley splitting, induced by the sharp 
interface. For large Un, a much larger gap emerges, A ~ 
200 meV, which cannot be explained by valley splitting, 
indicating an IS, as discussed below. 

Comparison between the three panels of Fig. [2] reveals 
that the TB model used in region III does not qualita- 
tively affect the nature of the IS. This is because the step 
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FIG. 2: (Color online) Ground state gap A as a function 
of the TB hopping parameters in region II. Results for the 
three barrier models (i)-(iii) are presented in panels (a)-(c) 
respectively. In each case, we assume an electric field F — 
0.05 V/nm, and a ID supercell that is large enough to avoid 
finite size effects. The circles in panel (b) indicate asymptotic 
points in region II, where the hopping parameters correspond 
to Si or Si02- The parameter a £ [0, 1] interpolates linearly 
between these points. 



in the CB is so large that the electron barely penetrates 
into region III. For this reason, and for simplicity, we 
adopt model (ii) throughout the remainder of the paper. 

Since the gap A does not exhibit strong features in 
Fig. [21 we consider two representative special cases, 
(u u = u S i, v u = Vsi) and (u u = u sio . 2 , v n = v S i0 2 ), as 
indicated by circles in Fig. HJb), and the line connecting 
these points. This line can be expressed parametrically 
through the equations un = (1 — a)ugi + ausi0 2 an d 
vu = (1 — ot)vsi + ctvsio 2 - We may then consider just 
one parameter, a £ [0, 1], which interpolates between a 
conventional CB state (a = 0) and an IS state (a — 1). 

Hybridization of IS and CB States. — We now demon- 
strate how IS can emerge and compete with valley states 
through hybridization, leading to large ground state gaps 
in the right-hand regions of Fig. [2j 

The five lowest energy eigenvalues of the ID TB hamil- 
tonian are shown as a function of the interpolation pa- 
rameter a in Fig. [3] [21]. On the left-hand side, we ob- 
serve nearly degenerate doublets. We refer to these as 
valley pairs because they share an envelope, but their 
fast oscillations are tt/2 out of phase. In this regime, it 
is appropriate to equate the ground state gap with the 
valley splitting Ay. Its value of ~ 2 meV is smaller than 
Fig. [2] because the electric field F — 0.01 V/nm is lower. 

On the right-hand side of Fig. [3[ an IS splits off from 
the CB into the band gap. In this regime, the ground 
state gap becomes quite large. However, the wavefunc- 
tions of the ground and excited states are very different, 
and they do not share an envelope, so it is not appropri- 
ate to speak of valley splitting. 




Interpolation Parameter a 

FIG. 3: (Color online) The five lowest energy levels in the ID 
TB model, for an electric field of 0.01 V/nm. The insets show 
the ground state electron probability densities corresponding 
to the representative values a = 0.05, 0.70, and 0.95. The 
region where strong hybridization occurs is highlighted. 



At intermediate values of a ~ 0.5 — 0.8 the IS hy- 
bridizes with the CB states, and we observe characteris- 
tic level anti-crossings, highlighted in Fig. [3] The mixing 
involves only the lowest state in each doublet, indicat- 
ing a valley selection rule. Such IS hybridization always 
enhances the ground state gap compared to Ay. 

The insets show the electronic probability density of 
the lowest ID eigenstate, obtained for three values of a. 
We observe characteristic CB and IS features, consistent 
with our previous discussion. The intermediate, hybrid 
wavefunction exhibits characteristics of both functions. 

Disordered Interfaces. — The hybridization of IS and 
CB states depends strongly on the unknown and locally 
varying hopping parameters of region II. The disorder 
can occur laterally, so that conditions favoring CB states 
or IS may coexist in a given heterostructure. Microscop- 
ically, the variations in hopping parameters can be inter- 
preted as the stretching and bending of Si-0 bonds to 
accommodate crystalline imperfections and the natural 
buckling of the interface. In some cases, an O atom may 
be completely absent, causing a missing link for charge 
transfer, which we model here as uu = 0. We now show 
that lateral disorder in the TB parameters may actually 
enhance the ground state gap, in contrast with interface 
height disorder, which suppresses the valley splitting [§]. 
We also investigate the possibility of lateral localization 
of a hybridized IS, in the presence of disorder. 

We introduce disorder into our 2D TB model by assign- 
ing random interpolation parameters a for each position 
y along the interface, in the range [—0.26, 5a], keeping 
5a constant in each realization. The lower bound of this 
range corresponds to the absent link condition, «n = 0. 
The upper bound 5a characterizes the degree of the dis- 
order. FigureQJa) shows the mean ground state gap A as 
a function of 5a. We conclude that any amount of bond 
disorder that can be modeled in this way will enhance A. 
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FIG. 4: (Color online) (a) The average ground state gap as a 
function of the disorder amplitude Sa, for 100 realizations of a 
supercell containing N y = 50 atomic layers in the y direction 
(2D calculation). The error bars correspond to a standard 
deviation of one sigma. The histogram inset shows the distri- 
bution of ground state gaps for the case Sa = 1. (b) The 2D 
electronic probability density |^(y,2)| 2 is shown for a "worst 
case" disorder realization. The dashed line indicates the in- 
terface position. All results assume F = 0.05 V/nm. 



When the interpolation parameter a ranges all the way 
to Sa ~ 1, the resulting gap is consistent with the large 
experimental value reported in Ref. 

The error bars shown in Fig.[5Ja) indicate the standard 
deviation in values of A for different disorder realizations. 
In the inset, we show a histogram for the frequency distri- 
bution of the ground state gap, for the disorder amplitude 
Sa = 1. In this case, all realizations give A > 15 meV, 
indicating that the bond disorder consistently enhances 
the ground state gap by a large amount. 

Large energy variations, caused by lateral disorder, 
potentially induce IS localization. Such localized states 
might not be visible in transport measurements like those 
of Ref. [5j. A detailed study of localization is beyond the 
scope of this work. However, we can check for lateral 
localization over smaller length scales. In Fig. HJb) we 
plot the 2D electronic probability density for a "worst 
case" disorder realization, where the ground state gap 
is particularly large, with A = 36 meV. We observe that 
the wavefunction spreads out uniformly across the lateral 
plane, due to its vertical extension far outside of region II 
where the disorder occurs. For a high quality Si/Si02 in- 
terface of width ~ 1 A this suggests that the vertical 
extension of the hybridized state may be able to overcome 
localization effects. This is in contrast with extrinsic IS, 
which are typically localized. Note that disorder in the 
lateral hopping parameters, not included here, may also 
contribute to localization. 

Conclusions. — We have shown that a typical ground 
state of an electron confined at a Si/SiC>2 interface is a 
hybrid of interface and CB states. The IS component 



causes a sizeable ground state gap, while the CB compo- 
nent is vertically extended, providing conduction charac- 
teristics similar to a pure CB state. Hybridization there- 
fore provides a plausible explanation for the measure- 
ments reported in Ref. 0, without invoking many-body 
physics. We also conclude that bond disorder may in- 
crease the ground state gap by enhancing hybridization. 
At Si/Si02 interfaces, this mechanism competes with the 
more widely studied suppression of valley splitting, which 
occurs at a rough interface The net enhancement or 
suppression of valley splitting then depends upon the de- 
tails of the heterostructure. 

Our results suggest that it may be possible to tune 
the ground state gap and the localization properties of a 
Si/SiC>2 interface by controlling the hybridization of IS 
and CB states in the ground state wavefunction. Con- 
trol parameters could include electric or magnetic fields, 
providing new directions for theoretical and experimen- 
tal investigations. More detailed studies of atomic scale 
Si/Si02 bond structures and realistic disorder models are 
also needed, to help pinpoint how the fabrication process 
may be used to engineer the ground state hybridization. 
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